Free expansion of a Bose-Einstein condensate in a ID optical lattice 



O. Morsch, M. Cristiani, J.H. Miiller, D. Ciampini, and E. Arimondo 
INFM, Dipartimento di Fisica E.Fermi, Universita di Pisa, Via Buonarroti 2, F56127 Pisa, Italy 

(February 6, 2008) 



c<i : 
o . 
o 

(N ■ 



O 



o 

> 
00 

in 
o 

(N 
O 

a 

B 

I 

O 

o 



X 



We have experimentally investigated the free expansion 
of a Bose-Einstein condensate in an array of two-dimensional 
traps created by a one-dimensional optical lattice. If the con- 
densate held in a magnetic trap is loaded adiabatically into 
the lattice, the increase in chemical potential due to the ad- 
ditional periodic potential is reflected in the expansion of the 
condensate after switching off the magnetic trap. We have 
calculated the chemical potential from measurements of the 
transverse expansion of the condensate as a function of the 
lattice parameters. 

PACS number(s): 03.75.Fi,32.80.Pj 



The properties of Bose-Einstein condensates (BECs) in 
lower-dimensional trapping potentials have recently at- 
tracted increasing interest. Both theoretical and exper- 
imental investigations have been aimed at studying the 
differences between ID |-§ and 2D traps 
with regard to the static and dynamic behaviour of BECs 
in such potentials. A thorough understanding of these 
properties is important as a prerequisite for predicting, 
e.g., the evolution of a BEC loaded into ID and 2D wave- 
guides. Condensates in 2D and ID have been realized in 
magnetic traps starting from a 3D situation by changing 
the aspect ratio of the trap and the number of atoms in 
the condensate Q|. In order to realize a ID condensate, 
optical dipole traps have been used to achieve the neces- 
sary asymmetry between the trapping frequencies [§,1). 
Similarly, 2D condensates can be created in an array of 
pancake-shaped traps provided by the periodic potential 
of a ID optical lattice 14|. For large lattice depths, 
i.e. in the tight-binding regime, the trapping frequencies 
along the lattice direction exceed those of the magnetic 
trap by orders of magnitude and the ID lattice hence rep- 
resents an array of 2D traps in which the motion along 
the direction of the array is frozen. Recently, Pedri et 
al. have calculated the variation in chemical poten- 
tial when the depth of the periodic potential is increased. 

In this paper, we report on experiments with magneti- 
cally trapped BECs loaded into a one-dimensional optical 
lattice and subsequently allowed to expand freely inside 
the lattice. After the creation of an array of 2D traps 
within the ID optical lattice, the magnetic trap is sud- 
denly switched off. The condensate is then free to expand 
in the radial directions whilst still being confined in the 
lattice direction. The measured expansion allowed us to 
infer the chemical potential and to test its dependence 
on the lattice parameters. 

Our experimental setup for creating Bose-Einstein 



condensates of rubidium atoms and loading these into 
ID optical lattices has been described in detail else- 
where Briefly, we obtain BECs of « 1 - 3 x 10'' 
rubidium atoms inside a triaxial TOP-trap which are 
then loaded into a ID optical lattice created by two inde- 
pendently controllable, linearly polarized Gaussian laser 
beams of wave- vector k and waist w — 1.8 mm, detuned 
by 25 — 40 GHz from the rubidium resonance line and 
propagating at an angle 9 with respect to each other. 
The lattice constant d of the resulting one-dimensional 
periodic potential is then given by d = 7r/(fc sin(0/2)), 
and in the following we shall express the lattice depth 
Uo in units of a re-scaled lattice recoil energy Erec — 
h^TT^ / {2m(P) corresponding to the respective lattice 
spacing d. In our setup we could realize both a horizon- 
tal lattice with d — 1.56 /im and a vertical lattice with 
d = 1.2 /im. In the vertical case, after switching off the 
magnetic trap we accelerated the lattice downwards with 
an acceleration a = 9.81ms~^ by chirping the frequency 
difference between the lattice beams in order to compen- 
sate gravity in the rest frame of the lattice. Apart from 
the different lattice spacing, the vertical and horizontal 
cases were, therefore, equivalent. 
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FIG. 1. Expanded condensate after 23 ms of time-of-flight 
without a lattice (a) and with a lattice of depth Uo — 20 Erec 
present during the expansion (b). The line plots represent 
the integrated profiles of the condensate density along the 
lattice direction; the reduced width p|| along that direction 
(horizontal in this figure) is clearly visible. The white ellipses 
indicate the initial size of the condensate in the trap. 



In order to load the condensate into the optical lattice, 
after adiabatically reducing the mean magnetic trap fre- 
quency uJtrap to the desired value we linearly increased 
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the lattice depth from to Uq in a time tramp = 150 ms. 
Since typicaUy the chemical potentials no of our conden- 
sates in traps with frequencies (Utrap/27r « 20 — 70 Hz lie 
between 50 Hz and 200 Hz, the adiabaticity condition fl^ ] 
tramp > h / fiQ was Satisfied even for small chemical poten- 
tials. After that, only the magnetic trap was switched off 
and the condensate was imaged after a variable time t^xp 
of free expansion inside the optical lattice. The effect of a 
deep lattice on the free expansion is clearly evident in the 
condensate images of Fig. ^ and in the measurements of 
the condensate dimensions along and perpendicular 
to the lattice direction {p±) shown in Fig. || denot- 
ing the half-width of a Gaussian fit to the density 
profile): Whilst in the lattice direction the condensate 
does not expand at all, the expansion in the direction 
perpendicular to the lattice is noticeably enhanced with 
respect to the expansion of the condensate in the ab- 
sence of the optical lattice. The lack of expansion in 
the lattice direction reflects the fact that the condensate 
has effectively been split up into several smaller conden- 
sates confined in the individual lattice wells, whereas the 
enhanced expansion in the perpendicular direction is ex- 
plained by the increase in the chemical potential when 
the lattice is ramped up. 

Taking into account the periodic potential, Pedri et 
al. have calculated that in the limit of large lattice 
depth [Uq ^ Erec) and assuming that the Thomas-Fermi 
approximation is always valid in the radial direction, the 
(local) chemical potential /ifc of well k is given by 

^Jik = \moJlS{kl-e), (1) 

where a;|| is the magnetic trap frequency in the lattice 
direction (similarly, lo±^ will indicate the trap frequency 
along the (observable) orthogonal direction). The well 
number k ranges from to the maximum /c,„ given by 

fc^^^ , ( 2^^trap \ ^'^ f 15 as d\ ^ 

where N is the total number of condensate atoms, aho = 
^Jfi/ {muj) the harmonic oscillator length, is the s-wave 
scattering length and a\\ is the width along the lattice 
direction of a Gaussian wavepacket in a single lattice well. 
We calculated a\\ by using a variational ansatz with a 
Gaussian wavefunction |Q . The total number of wells is 
given by n = 2fcm -I- 1. For the central well with fc = 0, 
the (local) chemical potential is 

^ , ^ 1/10 

where no is the chemical potential in the absence of the 
lattice. We have checked that considering only the cen- 
tral well for the perpendicular width reproduces to within 
a small error (< 1%) the result of a fit to the envelope of 
all lattice wells. 
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FIG. 2. Condensate dimensions (single shot values) p\\ (cir- 
cles) and (squares), versus the expansion time t^^^ in the 
presence of a horizontal lattice of depth C/o ~ 25 -Erec (filled 
symbols) and without lattice (open symbols). The conden- 
sate was initially released from a trap with ZZ7t,.ap/27r = 25 Hz 
(ljx/'^iI = l/\/2). The dashed and solid lines are the results 
of a numerical simulation (see text). 

We modeled the expansion of the condensate in the 
presence of the ID lattice using the equations derived 
in ||19|,po| with a slight modification along the lines of 
Ref. [[Zl: In the ansatz leading to the differential equa- 
tion of [|9| for the scaling factor A||(i — texp) — P\\{t — 
texp) I P\\{t — 0) along the lattice direction we replaced the 
atomic mass m by the effective mass m*(C/o), introduced 
in [Q in analogy with a solid state physics approach, as 
derived from a band structure calculation for a periodic 
potential of depth [/o- As can be seen in Figs. || and ||, 
taking into account the variation of p,k=o with Uq this 
model reproduces well our experimental data for the per- 
pendicular expansion of the condensate (the theoretical 
plots are corrected for the 5 /xm resolution of our imag- 
ing system). For the lattice direction, this approach gives 
the correct result for deep lattices (C/q ^ Erec) and re- 
produces reasonably well the qualitative behaviour in the 
intermediate regime (0 < C/q < 5 Erec)- We find, how- 
ever, that the experimental expansion along the lattice 
direction is considerably less than the theoretical pre- 
diction. We have checked that taking into account the 
finite momentum spread of the condensate when calcu- 
lating the effective mass only leads to a correction on 
the percent level and thus cannot explain the deviation 
of our experimental data from the numerical calculations 
neglecting mean field corrections. This might indicate 
that there are effects such as self-trapping [£1[ due to the 
mean-field interaction that further reduce the expansion 
in the lattice direction. 
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FIG. 3. Dependence on the lattice depth of the size of 
a condensate released from a trap with ZJt,.ap/27r = 25 Hz 
(LJx/i.tJ|| — \/2/l) after a time-of-flight of 22ms for a vertical 
lattice. For p± the solid line is the result of the numerical inte- 
gration of the scaling equations (see text) taking into account 
the variation of the chemical potential /ife^o with the lattice 
depth, whereas the dot-dashed line was calculated assuming a 
fixed po- In the lattice direction (circles), the experimentally 
measured expansion for p|| lies below the theoretical predic- 
tion. The error bars indicate the standard deviations of the 
averaged data points. 



For the perpendicular direction we have included in 
Fig. ^ the results of a numerical simulation without tak- 
ing into account the variation of the chemical potential 
with the lattice depth, i.e using fio instead of /ifc=o of 
Eqn. ^. Clearly, the resulting increase in perpendicular 
size is much less than what we observe experimentally 
and is related to the different re-distribution of mean- 
field energy when the degree of freedom in the lattice di- 
rection is frozen out. We have checked that the variation 
in lattice depth during the time-of-flight of up to 23 ms 
in the Gaussian beam profile (during which time the con- 
densate drops by up to 2.6 mm) does not affect the final 
widths. In fact, the mean- field explosion (during which 
mean-field energy is converted into kinetic energy) occurs 
during the first few milliseconds (equivalent to < 100 fim 
of free fall), after which the confinement along the lat- 
tice direction is almost perfect as long as the local lattice 
depth is larger than a few recoil energies. 

We could also directly deduce from our data the de- 
pendence of the chemical potential on Uq by calculating 
fJ-k=o from the initial perpendicular size p± (t = 0) of the 
condensate in the presence of the lattice inferred from 
the size p^^{t — t^xp) measured after an expansion time 
texp and the ratio X±{t = texp) / X±it = 0) of the scale 
factor X±{t). The chemical potential was then given by 



fJ-k=o = -muj^^p\{t = 0). 



(4) 



In Fig. ^ we plot the ratio iik=o /mo between the chemical 
potential as measured from the expansion of the conden- 
sate in the optical lattice and the chemical potential /ip 
in the magnetic trap alone. Agreement with Eqn 
very good. 
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FIG. 4. Variation with the lattice depth Uo of the chemical 
potential fik=o (in units of po in the absence of the lattice) 
calculated from the data of Fig. ^. The solid line is the theo- 
retical prediction of Pedri et al. 





40 -| 




35- 


E 






30- 


H 










25- 




20- 




p20 




-18 
-16 


JO 


- 14 


3, 


-12 




-10 





10 



15 
N 



20 



25 30x10 



3.6 
3.4 
07 3,2 
^3.0 
2.8 
2.6 



b) 




i 

H. 


■ 


- 


• 

• 


M 


H M 

W 

1 





10 



15 
N 



20 



25 30x10 



FIG. 5. Variation of the condensate widths (a) along 
(p||, circles) and perpendicular to the lattice direction {p±, 
squares), and of the the aspect ratio (b) with the number of 
atoms A^. The solid lines are predictions based on Eqn. ^ 
for the chemical potential /ifc=o in the presence of the lat- 
tice. In this experiment Uo/Erec = 20, iutrap/2n = 12 Hz 
with c^x/aj|| = V2/I. The error bars represent standard de- 
viations for the averaged data points at fixed final RF cut. 



We have also checked the dependence of fik=o on the 
number of atoms N in the condensate. As TV is varied by 
cutting into the (almost pure) condensate during the final 
stage of the RF evaporation, for a deep lattice the size 
of the condensate after a time-of-flight increases consid- 
erably in the perpendicular direction when N increases, 
as expected from the variation in chemical potential (see 
Fig. H (a)). In the lattice direction, however, the increase 
in size shows a weaker dependence on N than expected. 



For N > 2 X 10^, p|| even starts to decrease with in- 
creasing condensate number. Consequently, the aspect 
ratio p±/p\\ increases with N (see Fig. ^ rather than 
remaining constant as theoretically predicted. This be- 
haviour suggests that the condensate width in the lattice 
direction docs not fully reflect the change in chemical 
potential. When N gets large, mean-field effects might 
become important in the dynamics of the condensate ex- 
pansion, as mentioned above, and non-linear effects like 
self-trapping might reduce the observed width p\ | . 

Finally, we note here that in the limit of large lattice 
depths, our experiments effectively realize an adiabatic 
transformation between a 3D condensate and an array of 
2D condensates. The condition < TiuJiat of Ref. ^] 

(where ujiat = '^{Erec/K)\/Uf)l E^ec is the harmonic ap- 
proximation for the oscillation frequency in a lattice well) 
for the condensates in each well to be in the 2D limit is 
always satisfied for the small number of atoms in a sin- 
gle well (~ fO'^) present in our experiment. For an array 
of 2D condensates obtained by creating the condensate 
in the combined potential of the harmonic trap and the 
lattice, Burger et al. ||l^ have shown that in the case of 
their cigar-shaped condensate (with the long axis along 
the lattice direction), the transition temperature T^^ in 
the presence of the lattice is significantly lower than T^^ 
in the 3D case (i.e. in the magnetic trap without the lat- 
tice) . Calculating the critical temperature T^^ along the 
same lines for our system, we find that T^^ w T^^ due 
to the larger number of atoms per lattice site in our ge- 
ometry, and hence we expect no significant change in the 
condensate fraction in the presence of the lattice. In fact, 
experimentally we even find a consistently larger conden- 
sate fraction after ramping up the lattice. This result in- 
dicates that, with an appropriate choice of parameters, a 
ID optical lattice could be used to investigate adiabatic 
transformations between 3D and 2D condensates which 
could, e.g., be exploited to create condensates from ther- 
mal clouds by changing the dimensionality of the system, 
similarly to the change in the shape of the potentials in 
Refs 1^,^ for other geometries. 

In summary, we have studied the free expansion of 
a Bose-Einstein condensate in a one-dimensional opti- 
cal lattice. When increasing the lattice depth, an in- 
crease in the chemical potential of the condensate was 
measured, in agreement with a recent theoretical calcu- 
lation jlj]. The expansion along the lattice direction, 
though, seems to be strongly reduced in the regime of in- 
termediate lattice depth, necessitating a more thorough 
theoretical treatment taking into account mean-field ef- 
fects in the dynamics after the magnetic trap is switched 
off. Furthermore, the possibility of adiabatically chang- 
ing the dimensionality between a 3D trap and an array of 
2D traps should, with appropriate experimental param- 
eters, represent a way to reach quantum degeneracy by 
adding a periodic potential to a harmonic trap. The im- 
proved characterization of the behaviour of condensates 
in optical lattices is important with regard to adiabaticity 



criteria when a condensate is transferred into a periodic 
potential that could, e.g., be used to implement a neutral 
atom quantum computer. 
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